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Abstract 

We study the question of whether or not QCD predicts a pentaquark 
state Q + . We use the improved, fixed point lattice QCD action which 
has very little sensitivity to the lattice spacing and also allows us to reach 
light quark masses. The analysis was performed on a single volume of 
size (1.8 fm) 3 x 3.6 fm with lattice spacing of a — 0.102 fm. We use the 
correlation matrix method to identify the ground and excited states in the 
isospin 0, negative parity channel. In the quenched approximation where 
dynamical quark effects are omitted, we do not find any evidence for a 
pentaquark resonance in QCD. 



1 Present address: Department of Physics, Carnegie Mellon University, 5000 Forbes Ave.. 
Pittsburgh, PA 15213. 



1 Introduction 



The initial observation ^ of a baryonic resonance + with strangeness S = +1 
generated an enormous amount of experimental and theoretical activity. As the 
minimal quark content of the + is uudds, this was the first evidence for exotic 
hadronic states with more than three quarks, which have long been conjectured 
to exist |2] • There followed a number of experimental results ]3J consistent with 
a pentaquark resonance close to 1540 MeV, roughly 100 MeV above the Kaon- 
Nucleon threshold, with a surprisingly narrow width, possibly even as small as 
a few MeV . However, there have also been several null experimental results 
where the + resonance has not been found 0]. The uncertain experimental 
status of the pentaquark is an ideal situation for theoretical calculations to have 
an impact, particularly as the parity of the G + has not yet been measured (it 
is expected to have isospin 0). In fact, previous to the experimental discovery, 
it was actually predicted |2| in the chiral soliton model that a narrow resonance 
should exist at approximately 1530 MeV, which was completely consistent with 
the first observations. There have since been a number of attempts j?J] to model 
the behavior of the pentaquark, in particular its very narrow width, for example 
using quark models, Kaon-Nucleon bound states and QCD sum rules. 

The only known method to derive hadronic properties ab initio from a theory 
of strongly-interacting quarks and gluons is lattice QCD. Not surprisingly, there 
have been a number of lattice studies to see if the + can be predicted from 
QCD. The first attempts used standard hadron spectroscopy methods to explore 
pentaquark states in a number of isospin and parity channels, to determine for 
which quantum numbers could a pentaquark state be consistent with experiment 
[HI El IH9 ■ It was consistently found that the lightest isospin state is below the 
lightest isospin 1 state, but there was disagreement on the parity assignment of a 
possible + resonance. There were also studies which did not find a pentaquark 
resonance, only scattering states of weakly-interacting Kaons and Nucleons ^2 
112) . More recent studies which try to disentangle a possible pentaquark from 
Kaon-Nucleon states indicate that a pentaquark resonance is not seen ^3 ^] . 
However consensus has not yet been reached, as other work 1101115111?)] continues 
to favor a single-particle resonance consistent with the experimental + . The 
search for spin-3/2 pentaquarks on the lattice have also been conducted more 
recently by a couple of groups (E3E|); however, also without agreement on 
their conclusions. 

The purpose of this short paper is to try to add to the level of agreement 
among lattice QCD calculations. There are countless ways QCD can be reg- 
ulated using the lattice as a space-time cutoff, but all possible formulations 
should agree in their physical predictions when the lattice spacing vanishes and 
the continuum limit is reached. We use a particular improved lattice QCD 
action, designed so that physical continuum results can be extracted even on 
relatively coarse lattices. Because our action is also chirally symmetric at van- 
ishing quark mass (up to small parametrization errors), this allows us to reach 
light quark masses, which is not possible in many lattice formulations. We 
use the correlation matrix method to separate a possible pentaquark resonance 
from an interacting multiparticle system. As in the other studies, we use the 
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quenched approximation where dynamical quark effects are omitted, which is 
an enormous saving in computational difficulty. Although without justification, 
this approximation empirically gives qualitatively and even quantitatively the 
correct picture for the hadronic states whose width is narrow. The expectation 
is that this will also hold for any possible pentaquark states, if it is indeed as 
narrow as some experimental results claim. 

The paper is organised as follows: in Section 2 we give a short description 
of the lattice action and the details of the numerical simulation. In Section 3 
we discuss the correlation matrix analysis to separate ground and excited states 
in various channels. In Section 4 we present and discuss our numerical results, 
followed by our conclusions. 



2 Simulation details 

The lattice QCD action we use is based on the real space renormalization group 
|19j . The fixed point action is a close approximation to the exact renormalization 
group trajectory, so the theory has only a weak dependence on the lattice spacing 
and continuum results can be extracted even on relatively coarse lattices. If the 
continuum limit is reached earlier, this could bring a large reduction in the 
computational time required to work at a fixed physical volume. Fixed point 
actions have been constructed for 2- and 4-dimensional systems and in general 
have shown the expected good behavior |20| . In addition, the exact fixed point 
Dirac operator is a solution to the Ginsparg- Wilson relation [23, giving exact 
chiral symmetry at non-zero lattice spacing |22| . The fixed point action is one 
of the many possible lattice fermion formulations that have all the desirable 
chiral properties of the continuum theory, such as masslcss pions at vanishing 
quark mass. The action we use in the simulations is close to an exact Ginsparg- 
Wilson solution, hence the chiral properties are almost completely intact. 

The fixed point QCD action has been tested in a number of quenched QCD 
studies [21], including hadron spectroscopy, ir-ir scattering lengths, the chiral 
condensate, glueball masses, the topological susceptibility, unphysical zero-mode 
effects in quenched QCD, and properties of the Dirac operator eigenvalues. 
Large-scale hadron spectroscopy studies have shown the improved behavior of 
the action |251 126j . Hadron masses calculated at lattice spacings of 0.153 and 
0.102 fm have shown only a weak dependence on the lattice cutoff and are closer 
to continuum-extrapolated values than using the Wilson action at 0.05 fm or the 
improved staggered action at 0.09 fm. Pion masses as light as 220 McV could 
be reached without any indication of exceptional configurations (i.e. gauge field 
configurations for which the Dirac operator has unphysical accidentally small 
eigenvalues due to large fluctuations), which limit the lightest possible quark 
mass that can be reached. This is an indication of the good chiral behavior of 
the lattice action. 

In this paper we study an ensemble of 176 gauge configurations 2 of lattice size 

2 Four random configurations were irretrievable from the backup tapes making the total 
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18 3 x 36 which have already been generated for standard hadron spectroscopy 
|26| . From the Sommer parameter ro = 0.49 fm, the lattice spacing is found to be 
a(ro) = 0.102 fm, where the systematic error is at the % level and the statistical 
error is less than 1%. The gauge configurations are separated by 500 alternating 
Metropolis and pseudo-over-relaxation sweeps. The gauge configurations are 
fixed to the Coulomb gauge and Gaussian-smeared sources of width a /a ~ 2.3 
were used for the quark propagators. The ensemble covers a range of quarks 
masses, reaching pseudoscalar meson masses as light as raps — 300 MeV. 



3 Analysis 

The usual goal in standard hadron spectroscopy studies is to extract particle 
masses for the various possible quantum numbers. The typical method used is 
to calculate the correlation functions in Euclidean time of two- and three-quark 
operators with the desired quantum numbers. For example, the pion correlation 
function is CV(i) = (D,\0(t)0^ (0)\Q) , where the pion operator is O — V>75i/>. For 
large Euclidean time, the correlation function is dominated by the lowest energy 
state 

C(t) oc \(n\O\Q)\ 2 exp{-E t). (1) 

If the correlation function can be measured for sufficiently large time separa- 
tions, the higher energy states are suppressed and this is a safe method to extract 
the ground state energy. This can be seen in the effective mass 

m cS (t) = -\n[C(t + l)/C(t)}, (2) 

which develops a plateau when the higher states give only a negligible contribu- 
tion to the correlation function. However, if the energy separations are small, 
the higher energy states die out slowly and contaminate the lowest state. In 
addition, if the overlap of the operator with the ground state is accidentally 
small, it might be necessary to measure the correlation function for very large 
time separations to suppress the higher states. In these situations, the time 
separation required for the effective mass to develop a plateau might be unfea- 
sibly large. In principle, one can extract the energy of the excited states by 
including their contribution to the correlation function. However, in practice, 
this requires measuring C(t) to greater accuracy than is normally possible. 

A more reliable method to extract the energy states is to form a matrix 
of correlation functions of operators with the same desired quantum num- 
bers [271 |2Hj • For an n x n correlation matrix, the solution of the generalized 
eigenvalue problem 

Y.c^bf = \ n Y J c ij mf ) (3) 

between timeslices t = and t = 1 is used to construct the correlators 

C n {t) = Y J t )] Ci j {t)bf (4) 
i,3 

number of configurations 176 as opposed to the 180 which were used in Ref. 1261 , The conclu- 
sions are not affected by this loss in any way. 
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whose overlap with the n lowest energy states is optimal. The contamination of 
the ground state by higher states is minimized and the excited states themselves 
can also be extracted if the correlation matrix is measured to sufficient accuracy. 
The correlation matrix analysis is more reliable for extracting the ground state 
and it is necessary for getting higher energy states for multiparticle systems. 
If the + state does exist, its quantum numbers will be the same as a tower 
of Kaon-Nucleon interacting states. If the pentaquark is not the lightest state 
in that channel, at the very least the two lightest states must be extracted to 
positively identify the 9 + . 

To look for a possible pentaquark state + , some five-quark operator with 
the desired quantum numbers must be chosen. To construct a correlation ma- 
trix, more than one such operator is required. Without much guidance, it is 
difficult to tell which will have the best overlap with the + wave function. A 
number of operators have been tried, including ones based on a Kaon-Nucleon- 
type construction or using the diquark-diquark-antiquark structure proposed 
by Wilczek and Jaffe. Many possible five-quark operators can be related to one 
another by Fierz transformations and all of them must have some overlap with 
Kaon-Nucleon states. Even with the experience of the previous lattice QCD 
studies, it is not clear what the best choice is. To allow us to make a direct 
comparison with other work, we use the operator originally suggested in [J] and 
first used in a lattice simulation in [Sj 

9 = e abc [ulC l5 d b ]{u e s e i l5 d c T(u^ d)}, (5) 

as well as the color-rearranged combination 

KN = e abc [u T a C l5 d b ]{u c s e i l5 d e T (« <-► d)}, (6) 

where C = 7274 is the charge-conjugation matrix and the minus and plus signs 
correspond to the isospin and 1 combinations respectively. These operators 
have intrinsic negative parity, however, it is well-known that the correlation 
functions constructed from these operators receive contributions from negative 
and positive parity states. To project a particular parity state, an additional 
factor (1 ± 74) is included in the correlation function. 

The quark propagators were generated from Gaussian-smeared sources. Us- 
ing a smeared instead of point source has previously been shown to accelerate the 
decay of excited states in correlation functions and to generate effective masses 
whose plateaux extend to smaller time separation, allowing a more accurate 
determination of the effective mass. Smearing both the source and the sink in 
the correlation function has not in general shown much additional improvement. 
In our study, we use smeared sources and point sinks. This asymmetry means 
the correlation matrix is not Hermitian and so left and right eigenvectors must 
be determined in the generalized eigenvalue problem to construct the optimal 
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correlators: 



3 J 



^a^t a . (1) = A7i ^ a («)t c .. (0)) 



i i 



(7) 



4 Numerical results 

Previous lattice studies have examined a variety of isospin and parity channels 
to find a state compatible with the observed 9 + . With the exception of one 
study a possible signal of a genuine pentaquark state has been found only in 
the isospin negative parity channel. For the purpose of this short paper, we 
will concentrate on these quantum numbers. 

For this ensemble of gauge configurations, quark propagators were calculated 
for 10 quark masses, with pseudoscalar meson masses ranging from mps = 
300 MeV to 1390 MeV. Although one is ultimately interested in reaching the 
physical quark masses, it is also instructive to study the mass dependence of 
hadronic states. To calculate a correlation function of five-quark operators, the 
number of quark propagator contractions is much greater than for two- or three- 
quark operators. Hence we consider only some of the possible combinations of 
up-down and strange masses, m u d and m s , for this ensemble. In quenched 
QCD, where dynamical quark effects are omitted, gauge configurations with 
topological zero modes of the Dirac operator are not suppressed in the ensemble 
when the quark masses become small. The zero modes can give large unphysical 
C(1/to") contributions to correlation functions, which are only suppressed as 
the volume becomes large. This effect and possible solutions have been studied 
extensively for two- and three-quark correlation functions, but not yet for the 
pentaquark case. Since our physical volume is relatively small, roughly 1.8 fm, 
this might prevent us from extracting useful information from the lightest quark 
masses. 

We calculate the 2x2 correlation matrix of the operators Q and KN and 
construct the optimal correlators Cij (t) . Near the region where a plateau in the 
effective mass is observed, we fit the diagonal entries of the optimal correlators 
to a two-exponential ansatz 



giving us the ground and first excited energies Eq and E\. The best fit val- 
ues were obtained using the standard correlated x 2 method using the measured 
covariance to estimate the covariance matrix. The errors were obtained by a 
1024-point bootstrap sampling procedure. The different quark mass combina- 
tions were fitted independently and so correlations between the different quark 



Cu(t)=A 1 e- E ^(l + B 1 e- A ' t ) 



(8) 
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mass values were not taken into account. We should note that there are also 
backward-propagating contributions to the correlators, but these are suppressed 
by a factor of ~ 100 in the region where the two-exponential fits are performed. 
In Figures H and El we show typical effective mass plots and the correspond- 
ing two-exponential fits for the pentaquark ground and first excited state. A 
single-exponential fit can also be used, but only over a shorter time range, which 
limits the precision of the fitted parameters. The single- and two-exponential 
fitting procedures give completely consistent results, as shown in Figure To 
demonstrate the stability of the two-exponential fits, we show in Figure 01 the 
fitted energies for a particular quark mass combination as a function of the 
minimum timeslice included in the fits. The dotted lines are the % 2 per degree 
of freedom from the fits and the horizontal solid and dashed lines indicate the 
Kaon-Nucleon threshold and lowest non-zero momentum scattering states. The 
fitting procedure clearly works well. 

The meson and nucleon masses for degenerate m s and m u< i quarks have pre- 
viously been calculated for this ensemble and are tabulated in (25]. In Tables 
and [2] we give the results for the five-quark ground and first excited state en- 
ergies, and the pseudoscalar meson mass for non-degenerate quarks. Most of 
the data are for equal up-down and strange quark masses. The physical value 
of the strange quark mass is very close to am s = 0.078, so we also consider 
non-degenerate quark masses with the strange quark mass fixed at this value. 

Previous lattice QCD and experimental studies have shown that the Kaon- 
Nucleon interaction in the isospin s-wave channel is weak |29l \'M)\ . Therefore 
the energy of a Kaon-Nucleon scattering state with equal back-to-back momen- 
tum is closely approximated by 

E KN (p) = + \pf + yjM% + \ff, (9) 

where the momentum on a finite lattice is quantized p= 2irn/L. In Figurc|3wc 
plot the five-quark ground and first excited state energies and the three lowest 
Kaon-Nucleon scattering state energies as a function of the sum of the nucleon 
and pseudoscalar meson masses. Over a large range of quark masses, we see that 
the five-quark energies are in very good agreement with the weakly-interacting 
scattering states. As the quark masses are reduced, the five-quark excited state 
is more difficult to determine, but the ground state is still accurately measured. 
With m s close to the physical value and m u d < m s , the five-quark ground 
state agrees well with the Kaon-Nucleon ground state but we cannot make a 
strong statement about the excited state. The most important observation is 
that, for heavier quarks where our results are most accurate, we do not see 
any indication of a pentaquark resonance state, which should lie between the 
scattering state energy levels. The picture for light quarks in consistent with 
this, where the first excited state lies between the first and second excited Kaon- 
Nucleon scattering states, with much larger statistical errors. In principle, one 
would like to extrapolate the results to the physical quark masses. However, 
as we mentioned previously, quenched QCD correlation functions suffer from 
unphysical large contributions from topological zero modes at small quark mass. 
This complication might obstruct taking the chiral limit on this lattice volume. 
In addition, one would like to repeat this calculation for a range of lattice 
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Figure 1: The effective mass of the optimized pentaquark ground state cor- 
relator. The optimal effective mass as defined in Ref. 28 is also shown with 
diamond symbols. 



spacings to extrapolate to the continuum limit. However other hadronic states 
have shown only a very weak dependence on the lattice cutoff, so we expect this 
will not change our results significantly. 

Let us compare our results with other studies. The first attempts 3 [HI E] 
at pentaquark spectroscopy found an energy state slightly above the Kaon- 
Nucleon threshold in the isospin negative parity channel, which looked con- 
sistent with the experimental + . Further studies using Bayesian techniques to 
extract ground and excited states found energies consistent with K-N scatter- 
ing only Scattering states can also be distinguished from resonances via 
the volume-dependence of the energies, which can also be raised with hybrid 
boundary conditions ^2 • The amplitudes W n of the energy states contribut- 
ing to the correlation functions also have a distinctive volume-dependence for 
two-particle scattering. The correlation matrix method has also been used pre- 
viously to extract the pentaquark ground and excited state energies ITU ITfJ| . 

3 One of the referees of this paper has requested to emphasize at this point that the quark- 
exchange diagram between diquark pairs have been omitted in Ref. [gj. 
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Figure 2: The effective mass of the optimized pentaquark first excited state 
correlator. The optimal effective mass as defined in Ref. is also shown with 
diamond symbols. 



There is growing evidence that the isospin negative parity channel, the one 
thought most likely to describe the experimentally observed + , only has only 
weakly-interacting scattering states and no resonance. A stronger claim is made 
in ^2] i that a genuine pentaquark resonance should become a bound state at 
heavy quark mass, which is clearly not seen in any study. However, the picture 
is not completely consistent. In the pentaquark first excited energy state 
determined from the correlation matrix does not have the expected volume- 
dependence of two-particle scattering in small volumes. In [161 115j the ampli- 
tudes W n for the excited state appear to be independent of volume, signalling 
a resonance rather than a two-particle state, in contradiction to the findings 
in [U]. There is also the competing claim JOj that the + is actually in the 
isospin positive parity channel. 
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Figure 3: Single- and two-exponential fits to the ground state diagonal compo- 
nent of the optimized correlation matrix. For the single-exponential fit, the fit 
could only be done starting at large time separations. 



5 Conclusions 



Our conclusions are relatively straightforward. The correlation matrix technique 
is designed to extract a number of energy levels in a given channel. Previous 
evidence suggests that the physical + , if it exists, is in the isospin negative 
parity channel, above the Kaon-Nucleon threshold. In that case, the minimum 
requirement from a lattice study is to identify the two lowest energy states. 
Doing this, we can accurately determine the five-quark energies for a range of 
quark masses. We find that the ground and first excited states are completely 
consistent with weakly-interacting Kaon-Nucleon scattering states, with no in- 
dication of a pentaquark resonance which should lie between them. Although 
we work at a finite lattice spacing, previous work has shown that this lattice 
action reproduces extremely well continuum properties of hadronic states, so we 
do not expect our observation to change significantly in the continuum limit. 

As the title of this paper suggests, absence of evidence is not evidence of 
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Figure 4: The fitted energies for the pentaquark ground and first excited states 
and the quality of the fits as a function of t m i n , the smallest time separation 
included in the fit range. The horizontal solid lines are the Kaon-Nucleon thresh- 
old (including the error) and the horizontal dashed lines are the energy levels 
of Kaon-Nucleon scattering states. 



absence and there are many caveats one can include. This and all other lattice 
QCD pentaquark studies have ignored the effect of dynamical fermions, which is 
an unjustified approximation. The most accurate results ruling out a pentaquark 
resonance are at unphysically heavy quark masses, requiring an extrapolation to 
the physical values. We work at a finite lattice spacing and do not extrapolate to 
the continuum. However we do not expect these effects to qualitatively change 
our conclusions. We believe a more serious issue is our choice of operators in the 
correlation matrix, which may have little overlap with a genuine + state, if it 
exists. We cannot rule out this possibility, and there may even be other reasons 
why we do not see a pentaquark resonance. However we can only conclude that 
we find no evidence that QCD predicts the + state. 
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Figure 5: The fitted pentaquark ground and first excited state energies as a 
function of the sum of the nucleon and pseudoscalar meson mass. The degener- 
ate quark mass data are the dark blue and red circles, the non-degenerate quark 
mass data the light blue and orange squares. The hollow symbols are the three 
lowest Kaon-Nucleon scattering state energies. The dashed lines indicate the 
physical values of Mr + Mm- 
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Table 1: The fitted energy values for the ground and first excited states. The 
two states were fitted separately to two exponentials after the diagonalization 
of the 2x2 correlation matrix. 
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0.2817(15) 


0.037 


0.078 


6 


11 


1.03 


0.39 


0.2884(17) 


0.045 


0.078 


5 


12 


1.41 


0.21 


0.3031(14) 


0.100 


0.140 


5 


10 


1.16 


0.33 


0.4656(15) 



Table 2: The pseudoscalar meson mass for non-degenerate quark masses. 
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